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ABSTRACT 


The difficulties encountered when implementating ortho- 
graphic reprojection of 3-dimensional image data Omemuee 
2-dimensional screen are considerable. They arise princi- 
pally because of the size of data being manipulated and stme 
tendency for underlying or overlying structures to obscure a 
clear view of the desired image. In this work, implementa- 
tion was performed with an orthographic reprojection tech- 
nigue and many heuristic approaches were used to resolve 
some of the related problems. 

Several coordinate rotating algorithms were tested in 
this work. Among them the Precalculation and Indexing 
method proved to be the most efficient algorithm. 

Due to the disparity of the viewing-coordinate grids and 
the voxels of the volume data after rotation, 3-dimensional 
interpolation is required for appling the reprojection tech- 
nique. Several methods implementing linear interpolation 
have been tried. Interpolation with a cone-shape kernel is 
the most appropriate method in the 2D situations and can be 
easily extended to a sphere-shape kernel in 3D Situations. 

The orthographic reprojection method includes a single 
plane dissection capability. Results on an artificial tesge 
data are collected using the above algorithms. 

Several resulting images are included as well as’ the 
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I. INTRODUCTION 


A. GENERAL 


The representation of 3-dimensional(3D) Sbyects on 2 
2-dimensional(2D) screen presents an interesting problem in 
computer Science and computer engineering. The application 
meee this technology is common in areas such as computer 
graphics, CAD/CAE, and medical imaging. The study here will 
be concerned with amore specialized area of showing 3D 
objects on 2-dimensional (2D) Screen. A variety of 
techniques have been developed for this purpoSe. 

A discrete set of 3D volume data consists of a collec- 
tion of rectangular parallelepipeds. These are referred to 
as volume elements, or more commonly voxels. Each voxel has 
an associated value which represents the average intensity 
of a parallelepiped referred to as the density value. The 
density value is an integer from a finite set on. Dee). 
If Dyin =9 and Dmax=l, we say that the image is a binary 
image. ) 

Pemercetiisertil “to Classify 3D display techniques into 
different categories. [Ref. 1]. The main criterion used to 
classify them is whether the technique displays a whole 
scene or a only cross-section of a selected plane in the 3D 
volume data. 

WHOLE SCENE DISPLAY TECHNIQUE---The entire scene of 3D 
data can be displayed by the use of a vibrating mirror. A 
Single slice of a 3D object can be displayed onto a 2D 
Screen at a short time. If we then display many such slices 
in a rapid succession and view the screen through a 
vibrating mirror of variable focal length, we can make the 
Slices appear to be in BENGE corre ck slopene cil) ie yeksan erg’ 


jmemative to each other. 


Ieae 


Another way of displaying the whole scene is to project 
it onto the 2D screen. The work studied here is related to 
Chis@4approgden. 

SURFACE DISPLAY TECHNIQUE---Using the surface  disieea 
technique, we are interesting in dipicting the appearance of 
surfaces of an object in the scene. This involves 
sepment ing the original volume data into an object anda 
background, for example, by creating a binary image by which 
a lis asSigned to a voxel in the object and a 0 to aver 
in the background. Various techniques such as 1D Gomme 
based, 2D surface based, and 3D parallelpiped based or 
Sphere based display techniques are examples of surface 
displays. 

REPROJECTION TECHNIQUE---There are two kind of teehee 
niques in the reprojection method: Orthographic reprojection 
and Radial reprojectiton., (Rete. Orthographic reprojec- 
tion is performed numerically in the computer by summing the 
values Of voxels along parallel paths through the 
reconstructed computed tomography(CT) scan image. 

Radial reprojection is an alternate scheme in which the 
voxels of the CT Scan reconstructed . volume are projected 
onto a cylindrical surface surrounding the object, "Pamives 
than orthegraphically onto a plame- 

The reprojection of volume image data onto a 2D Screen 
sometimes cannot provide sufficient information because the 
size of any given structures 1S not dependent on the 
distance from the observer. Perspective does not exist in 
the orthographic reprojection. Therefore it is necessary to 
rotate the volume image data to the desired viewing-angle 
before reprojection. In addition to the rotation, numeiiee 
dissection and numerical dissolution techniquesma- intro- 


duced to show selected portions of the reconstructed volume. 
PRete. 2 |. 
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THE SCOPE OF THE STUDY---The study of techniques of 
Babin a ine eEepLojeectiom et 3D objects after rotation to a 
Polceredmev1ewine-anele is= the focus in this work. The 
computer technique for single plane dissection is also 
incorporated here. Generally 3D volume data consists of a 
huge amount of voxels. Therefore, an efficient algorithm is 
essential for the manipulation of volume data. 

After rotation to the selected viewing-angle, each voxel 
of the volume data does not lie on the grid of the viewing- 
coordinate system. In other words, each voxel looks like a 
polygon rather than a parallelepiped along the viewing-angle 
from the display screen. This will yield an obScuring effect 
between voxels during projection. Therefore, interpolation 
is necessary to resolve the contribution of voxels to each 
pixel in the viewing-plane. 

tee miaitieeempnasis im this work 1S on the study of 
efficient rotation algorithms and the interpolation method. 

ime Chapter Ie various EOtatvivonu a leorttnms are 
discussed, as well as the results obtained through the use 
Meeartifical test data. 

In Chapter III different interpolation methods using 


various kernel functions are presented. 


ie LeCHNIQUES OF ORTHOGRAPHIC REPROJECTION 


1. Reprojection Image Generation 


[Rew precess ~of reprojection tome lust raved in 
Figurel.1l. A linear array of squares on the viewing-plane 
corresponds to the picture elements(pixels) at the level k 
of this cross-section. The intensity value of the pixels in 
the array is the sum of those voxels along the projection 
path which are perpendicular to the projection plane. When 
all volume elements at the level k have been projected then 


the voxels at level k+l are projected. 
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Potation of Volume 


Viewing-Plane 


Reconstructed 
Voxels at 
Level K 





Reprojected Pixels 
at Level K 


Four emer Reprojectvom Processes 


The resulting 2D array of pixel values Gana 
rescaled and displayed onto the specific display device. 


Subsequently it is viewed as an image on a TV monitor. 
2. Numerical Dissection 


The dissection technique is used to reduce the 
effect of obscuration and to clearly show an internal Stige. 
ture of the volume data. In this technique, a plamemiou 
dissection 1s highlighted by selectively dimming the voxels 
in front of the plane and ignoring all the voxels behind the 
plane. Dimming is accomplished by replacing the value of 
each voxel with the product of an original voxel dense 
value andaeconstant less than unity. For example, a 
constant 0.1 would result in a intensity reduction tomes 
of its original value. | 

The result is that the internal structures "Jae 
plane of dissection are more clearly visible. It  i2S5aieaa 
easy to see the spatial relationship of the object plane to 


those structures in front of the plane. 
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See Gmetted lr yis so luUELOn 


Selective numerical dissolution iS a process whereby 
bme relative contributions of selected voxels can be 
decreased before rotation. fier samanner obscurine st ruc- 
tures are only partially removed and this results in a 'see 
through’ effect. 

The criterion for numerical dissolution is based on 
the density difference between the structures within the 
reconstructed volume. The density threshold is identified 
empirically by first Choosing an intitial threshold and then 
adjusting the value to achieve the desired dissolution 
effect. Once the deSired threshold is identified, numerical 
dissolution is accomplished by dimming all voxels with 
values below threshold during reprojection. The dimming 


process is the same as that in the dissection process. 


PS) 


IIT. COORDINATE TRANSFORMATION 


A. GENERAL 
1. Int rodgetion 


Several methods have been proposed for the display 
of 3-dimensional(3D) information contained within a4 ste 
of parallel Computed Tomography(CT) cross sections. ))iimemeee 
referred to as a 3D or volume reconstruction problem. “The 
reconstructed volume data Waele be displayed on a 
2-dimensional(2D) TV screen uSing reprojection, dissection, 
and dissolution “rechniaues: Although the problem of obscu- 
ration exists in the reprojected image, it iS not as Severe 
as inthe case of conventional radiographs because the 
reconstructed volume can be viewed from different viewing- 
angles. Before reprojection the object can be mathematically 
oriented to a desired viewing-angle. The choice of desired 
viewing-angles will enable radiologists to discover the 
important information from the image. 

The objective of this chapter 1S to present an 
example of 3D°Goordinate transtormacaon: Because digital 
image data is comprised of a huge amount of 3D data points, 
our rotation algorithm must be efficient. AlSo, cCoOoTrd#iiaae 
rotation invokes time-consuming floating-point operations 
the computers. An inefficient algorithm will prove unwieldy 
and lead to an undesirably long processing delay before 
results are available. In the final part of this chapters 
hardware alternatives to the software design which jaa 
implement the rotation algorithm developed in this chapter 


will be considered. 
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Zaeecenter of Rotation 





MicMeecnrCummememrTonatlOn issa@syoitrary, although it is 
generally chosen to be near the center of the object (Figure 
ZeMeniniusmicmcone So that after rotation the reconstructed 
miaee willesnot be shrrted too much out of the effective 


viewing-window defined on the screen. 
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(A) (B ) 


Figure 2.1] Center of Rotation(A) and Viewing-Plane(B). 


Pee Robetiom Types 


farea kimads Gt 8rotation are used for this work: 
Rotation about X-axis (Elevational), PORAETON about Y axis 
[Merizomtal), and rotation about both-axes(a combination of 
rotation about the X-axis followed by rotation about’ the 
Y-axis). These three kinds of rotation will be enough to 
provide a desired viewing-angle from any directions. Figure 
2.2 illustrates examples of elevational and horizontal rota- 
tions. Figure 2.2(A) shows +90 degrees rotation about X-axis 


and Figure 2.2(B) shows +90 degrees rotation about Y-axis. 


LZ 





(A)90 Rotation (B)Original Image (C)90 Roetatran 
SOU t + Aoi about Y-Axis 


Figure 2.2 Direction of Rotation: 


Rotation about both-axes is performed by the combination of 


elevational and horizontal rotations sequentially. 
4. Coordinate System 


Figure 2.3 shows the coordinate system we will use 
Ehroughout- this. works 
A coordinate system is selected according to the right 
hand rule. As shown in Figure 2.3, the horizontal line 
represents the X-axis, the vertical line represents’ the 
Yass and depth is labelled Waite i the Z-axise 
Object-coordinates are those we define relative to the 
object during image sampling, and viewing-coordinates are 
those defined for displaying the image. The viewing- 
coordinate system always corresponds to the coordinadtemimm 
display screen. Thus the orientation of the object with 
respect to the viewing-coordinates defines the rotate 


angle. 
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es >Object-Coordinate 


xX, 
Xe Vae2 Vrewine -woordinate 





Figure 2.3 Coordinate System. 


B. COORDINATE TRANSFORMATION 


ie oo hbip li cation 


Usually weacoondmnate  transitormation 1s done by 
[meciolying the original coordinate by a transformation 
matrix. DemvemUSec matrix motation this will result in a 
meetor transformation from an R*® space to another R*® space. 
We can define a function T: R? --> R*® by T(X) = AX where A 
Mmema transtormation matrix and X iS a input vector. Observe 
meveett Seis a 3xl matrix and A is 3x3 matrix, then the 
Peeocuct AX is @lso a 3xl matrix. hhas — maese. wage tio: R* . 
Mis kind of linear transformation is called ‘matrix trans- 
formation . melcehpratediiy aes) to 3D Eransftormation matrix 
meamires 3 basis functions. Suppose we want to rotate a 3D 
vector X = xx + yf + 2&8. Mics transtOrmation matrix for 


fers Operation should be a 3x3 matrix. In matrix notation we 
write 


ipa) Dube x | x= Agim y (dt eZ 
y dame fii drt rey + £2 
Z Sa me et Zz Pet My + 12 
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If the determinant of A is 1, it produces 4a pure rotageem 
about the origin. Before considering rotation about aneeeeue 
trary axis, two special cases are examined: rotation about 
the X-axis (elevational ) and rotdevor about the 


Y-4xe (hortzenealee 
2. Rotation about %-2xi06 


In this case the rotation matrix will have @zeaeaee 
the first row and first column except a unit value aon the 
main diagonal. The other terms are rotation dependent and 
will be determined by considering the specific rotation 
angle of a point. Figure 2.4 shows us this situation geome- 
trically. Rotation about only one axis is the Same "sua 
transformation because the rotational axis is fixed in 
Space. Here the rotation angle is defined to be positive in 


the counter-clockwise direction when a vector rotates 


= 


~~ 
= In this £1e¢ureser 1S ane uneve ctor 





Figure 2.4 Rotation abowt X-Axwee 


Im the Pie 2:4, thesx eee and Z represent the object- 
coordinates and X*, Ye, and Z* the viewing-coordinate system 


respectively. Thus our matrix transformation fomula is 
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Ge, y ye) = ea Eaenee sy 2 iz ° 


xe) = fl 0 0 x = x 
y° 0 cos@ =sin6| ly ycoos§ - zsin§ 
Ze 0 sing cos§ Zz VYSINGemezcos OC 


In figure 2.4, we let x = x, y = cos#, z= sing. After 


rotation we have 


=) = xX , 
ye = cos(§+ od) 
Ze = sin(@+ ¢) 


cosPcos®@- sinOsing= ycos9 - zsin§@, 


Moco Meameasosimg, — ysing + zeos@. 


iMiiomEnestiitiumeas Just the “matrix multiplication, the 
maaeaiet Of a transformation matrix and a vector in the 


object-coordinate. 
Dee Ret abLOnN about Y-Ax1S 


Rotation about the y-axis can be shown in the same 


Manner as rotation about the X-axis above. (See Figure 2.5) 


POS 5 Ae = IEC RO 520 

xe = cos § 0 sin6] | x = xcos@ + zsin§ 
y° O il O y y 

z° -sing 0 eos ez -xsinQ + zcos@ 


4. Rotation about Both-Axis 


Rotation about an arbitrary axis requireS a more 
complex operation. But we can restrict ourselves’ to two 
rotation directions for adequate display purposes. Ginaes 


MMm~olnaclon Of two directional! rotations can be considered 


fmqeecascaded Operation, first about the X-axis and then 
about the Y-axis or vice versa. Pie aocat rom about © the 
heaxis, then rotation about Y-axis 1s chosen, then the 


Mmereatlon matrix will be 
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Figure 2.5 Rotation about Y-Axise 


i 0 0) e0sr Of sim: = COSsa O sinr 
0 cos -sing Oo 1 0 sin@sinr cos -sin§fcessr 
O sin cos§6}ltsinr O cosr -cosQ@sinr sin§ cos§cosr 


This matrix could not be applied in the reverse direc- 
tional rotation, because a matrix multiplication bee 
generally not equal to BA. Therefore we have to be careful 


in our selection of matrix multiplication Sequences. 


C. PROGRAMMING METHODS 


General procedures for the image display are as follows: 
ee Set up the object-coordinate value, (2) Caleulatemiane 
corresponding viewing-coordinate value by coordinate trans- 
formation, (3) Access a pixel and place it in the Gauges 
Sponding viewing-coordinate, (4) Perform projection @eame 
Gas pila 

The following algorithms use different methods for coor- 
dinate transformation. Because coordinate transformation is 
the central issue in this problem the program efficiency 


will mainly depend on it. 


af 


Plies Cima c ree Ltiplicateron 


Suppose we want to rotate a 3D object. Every data 
point in the object-coordinate has x, y, z coordinates which 
form a 3D position vector. To rotate this object we have to 
multiply the position vector by a 3x3 transformation matrix. 
As shown in Table 1, the transform of each vector requires 9 
real multiplications and 6 real additions. If one dimension 
of the volume is N, the total calculation reaches 9N? multi- 
feesecations and ~6N° additions. ities tkow senart for imple- 
menting this method is shown in Figure 2.6. In a PASCAL 
environment, one record consists one plane and each record 
is accessed sequentially. Whenever an image plane is 
PeccsccGeeene | propram notates the coordinates of pixels 
contained in the plane one by one. Therefore the subroutine 


Mimieatioain § 1s called N* times. 


Define Viewing-Plane 
Open Data File 


Get the Direction of 
Rotation and the Angle 
Calculate Elements of 
Transtormation Matrix 


Read One Image Plane 
yeeliults |For Each Pixel, Transform a Object- Ve 
6 Adds {Coordinate into a Viewing-coordinate , 
Graces 
Times 


Value 


lBie 
N 
a 
iterates 
N* Times 
Close-Data File 


Display 


Paeure 2.6 Miewciabe wot the Direct Matrix Multiplication. 
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2. Decomposition with Fixed Coerdimage 


Let uS conSider a rotation about only senewacne AS 
shown in Equation 2.1 the transformation matrix ine iia 
zeros and aunit value in the main diagonal so the center 


axis of the rotation wilPwniot bDeeehameq 


xe) = (1 0 OY(x) = + y-O0 + 2-0 

y° O cos®8 -sin®O/ly x-0 + yrcos@ - zesin§] (2.1) 
z° O sinQ cos@J[z x°0 + yesin§ + zecos8€ 

sea) = Xx 

ye ycos® - zsin§ (2ae2,) 

ze zcosQ@ + zcos§G 


Equation 2.1 shows an example of rotation ab@u@uaee 


X-axis uSing the direct matrix multiplication Mme@ieus 
Equation 2.2 Suggests a more efficient implementation 
method. If the matrix multiplication is simplifleqd@tgum 


equations as shown in Equation 2.2, it resultS in an ecConame 
of 94 mMultiplzeations and 4 additions for each pixem 
transformation. 

Another point in comparing the two equations 1S )eiae 
the first row of the transformation matrix conSiStS Giese 
a unit value and two zeros. ThisS row maps a incoming vector 
into itself. (This is the fixed coordinate.) The other two 
coordinate values will be changed according to Equations 
This method is represented in the flow chart shown in Figure 
Zed Oe The inner loop does not include a call for a tai 
multiplication operation, but the outer two loops cai==eme 
Subroutine ‘Rotation’ which includes matrix multiplvea@euem 
Operations. Therefore the total number of calculatalene 
required are 4N* multiplications and 2N? additions. 

This method of decomposition with fixed coordinates 
is not very helpful for the case of arbitrary direction 


Ob ae 1 Gin. The transformation Sia bias of the arbitra 
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Chijeeectlomal rotatiecnmemas Yonly ene zero in the first row so 


the economy of calculations will be very small. 


Define Viewing-plane 
Open Data-File . 


Cetmune Direction ot Rotation and 
the Angle 
Caen lavewsilements oF 
Tans Pemmaiad Oflema tid x 


Read one Image Plane 


G4 Mults | Transform Each Imcoming Pixel (" 









2 Adds Lig Ete) Nene We Bacon ae Value By (‘ Iterate 
Seas ee ess ecomposed Equations Times 
ot Two semaneine coordinates pec rat 

imes 


AsSign Fixed-Coordinate Iterate| 
N Times 
Close Data-File 
Paisp lay 


Figure 2./ Flow-Chart of Decomposition with Fixed 
Coordinate(Rotation about One-axis). 


oe Precalculation and Indexing Method 


Another efficient method makes use of indexing 
instead of multiplication. This method is represented in the 


following equation 


faap= fell cl2 e113 a el ate? y+ Gl 37 
ye eZiege22 «e623 y ene? yout C2 37 (243)) 
ze eoamej2 6633 Z. Cia mes ly + C332 
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The elements (Cj) of the transformation matrix are fie 
after the rotation angle 1S specified, and the x, yj) eee 
values of the object-coordinate iterate individually. We 
can calculate the multiplication of each Coordinate wae 


corresponding matrix element ( Cuxk) before accessing pixels 


of volume and then store those values in an array. After 
this, by indexing the array according to the incoming 
object-coordinate and adding these components together, we 
can calculate the viewing-coordinate Valier Implementing 


this idea can allow us to avoid the time-consuming floating- 
DOPhit Mult lp wheat Toms. It 1S represented in the flowchart 
of Pupure 2.6. As shown in the figure the Indexing method 
does not use a ‘Rotation’ subroutine. 

The array contents also can be calculated by succes- 
sive additions as shown in the program ‘Rotation_3 Dim’. 
When the multiplication of a coordinate and matrix-element 
is calculated, a coordinate value always increases by one 
Write. For example, XK: = XK + 1, so the multiplicatwem 
(CreXy ) becomes Cy X pay = CX et Gy, CyXpe2= CyXky + Cysseretc. 
Only the first term requires multiplication of 4 pairs of 
elements. The total number of additions requi1Tredsiswyas 
+2N7. The 4N additions are the effective multiplicatiomomem 
the objective coordinate and transformation matrix element. 
The remaining 2N* additions result from calculations Gimeae 
viewing-coordinate value by means of double iteration loops. 

Because all floating-point multiplicae tomes are 
avoided by indexing and adding, this method results in large 
processing time Savings. The total number of calculations 
required is shown in Table 1. Compared to the other two 
methods, Table 1 indicates that this method is the most 
efficient of the three. But this method requires some 
working space in memory for storing the precalculated 
valves. In a VAX-11/750 floating-point notation (where a 
floating point value occupies 4b rece 12N bytes are 
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Becuiemed@Lor a  Gotatbion about ome axis, ~and 24N bytes are 
required for rotation about both axes. This is a reasonable 
memory requirement considering the processing time Savings 


oo 


realized by using this method. 


Define Viewing-Plane 
Open Data-File 


Cop aim treet on UrLenotation 
and the Angle 


SeMultes \Caleulare ene Mulltiplications of 
4N Adds | Coordinate and Matrix-Element, 
and Store Them in an Array 


Read One Image Plane 


Romenacn Pixel, Transtorm the 
Object-coordinate into the 

lewing-Coordinate Value 
by Indexing the Array 






lterate 
Assign Fixed-Coordinate Value ( N Times 


: (ee iealice 
Close Data-File 
Duss p lay. 


Figure 2.8 Flow Chart of Precalculation.and Indexing 
Method(Rotation about one-axis). 


Pee EROGRAM IMPLEMENTATION 


The program PRoie ss Dam which implements the 
Peeealculation and Indexing method is included in Appendix 
A. In this PASCAL program, ‘Objectcd' and 'Viewcd' repre- 
sent the object-coordinate and the viewing-coordinate 


respectively. Data tcc wro Nirccl the program oa Tl 
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artificially created double pyramid( N = 64 ). This 3Digeeee 
is shown in Figure 2.9 which portrays the Shape of ape. 


connected double pyramid. The density of the test data 





Figure 2.9 Double Pyramid. 


Sane 


- 
a gata { 


1 cocoa i 
Viewlng-444 | 512 
joe femme 


Full COMTAL Screen 





(A) (B ) 
Figure 222 00 Test Data(A) and Display Butierteee 
increases toward the center of pyramid. The object coordi- 


nate is defined so that the center of object is located at 


the center of the coordinate. The size along one dimension 
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of the data is 64, but the program needs 111xlll elements 
Size of a viewing-plane. In the worst case one side of the 
reconstructed image can have’ the diagonal size of the 


original image(Figure 2.10). 


E. NUMBER OF CALCULATIONS AND PROCESSING TIME 


The program ‘Rot 3 Dim' which manipulates our test data 
memeertin On the DEC VAX-11/750 computer. Results are 
displayed on the COMTAL image processing system which itself 
has a PDP-11/23 processor. Table 1 shows a comparison of 


program efficiency and execution time among these 3 methods. 


TABLE I 
foie bee Comparision of Methods 


Rotation about Oise hho mera otlt 
One-Axis Both-Axis 
Berect Matrix |9N* Mults over 25 9N? Mults jover 10 
Mettiplication|6N* Adds secs 6N? Adds mins 
Decomposition |4N? Mults 























with Fixed 
Coordinate aN “me Adds 













Precalculation|/4 Mults 

and Indexing |4N+2N? Adds 
(16N Bytes 
Storage) 






“ieee bine) sazerot dimension. 


Pees ceaie varue calculation. 


F. HARDWARE REALIZATION 


Modern computers have only one CPU. In this case each 


line of code must be executed sequentially. But digital 


Ze 


image data 1s composed of a huge number of pixels. Therefore 
a special purpose computer having parallel processing capa- 
bility is becoming popular. Here two simple hardware designs 


are suggested to realize the above algorithms of coordinate 


rotation. This hardware decreases the burden of iterative 
Matrix mle tpleve ad elem in the coordinate rotation Wagga 
releases the CPU to do more importame baci Also 1t Wag 
increase the speed of the image rotation These two design 
schemes use different components. Figure 2.11 illusteuape. 
the first hardware design scheme. The first design consists 
of 3 adders and 9 floating-point multipliers. An in@Omaimeee 


coordinate is multiplied with the first column of the tragee 
formation matrix, a Y coordinate the second column , ance 
coordinate the third column respectively. Then the multi- 
plied coordinates and matrix elements are added to make 
viewing-coordinate values as shown in Figure 2.11. 

Figure 2.12 shows the second hardware design. The second 
one consists of 3 adders and 9 / 4N bytes memories instead 
of multipliers. Each set of the monica has the capacity of 
4N bytes and individual memory element has 4 bytes width. 
The precalculated e(staneeelisn ieee and matrix elements will be 
stored in the memories. These components are indexed 
according to the incoming object-coordinate values, and then 
added to calculate viewling-coordinate values. It is clear 
that an indexing is faster than a floating point multipass 
t20n, As may be expected there is a trade-off between cost 


and speed. 
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Memory ui 
Goya } 


Bieure 2.1) Pikestemocmeme( 9 1 shlodtime- Point Multipliers ). 


Balt 





Bipure 2. 12 Second Scheme( 9 / 4N Bytes Memories ). 
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ition INTERPOLATION 


A. LNTERPOLATION IN IMAGE PROCESSING 


Interpolation is the process of estimating the interme- 
diate values of a continuous Signal from discrete samples. 
Interpolating 1s used extensively im Wdweea ta.) image 
processing to magnify or reduce images and to correct a 
spatial distortion. ery 3-dimensional (3D) digital image 
processing, each volume element(voxel) value represents the 
intensity value of a rectangular parallelepiped. eta S 
rectangular parallelepiped is not necessarily a cube, but it 


1s assumed to be a cube that has equal length edges for our 


purposes here. A digital intensity value 1S generated by 
mrarhoe-to-dipital conversion of a continuous intensity 
Ssienal. Therefore, a digital data produced in this way can 


involve a certain amount of aliaSing or blurring in the 
digitization process. But this effect is so small that we 
Sam ienore it. An exact image restoration will thus be 
possible if an appropriate interpolation method is applied. 
The goal of the interpolation of an encoded medical image is 
to reconstruct the original image as closely as possible. 
Because of the great amount of data associated with a 
digital image, an efficient interpolation algorithm is 


essential. 


jae INTERPOLATION METHOD 
ee One-Dimensional Interpolation 


An interpolation kernel function is a special type 
@eeapproximating function. A fundamental property of an 


maitetpolation kernel function is that it must have the 
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S(X) = SQUARE ****** square, 9 =) comvonieamen 
iene s 


-_— aap ow oe Go oo GD 


(A)Nearest Neighboring(n=1) (B)Linear Interpolation(n=2) 


(C) Cubic Spline(n = 3) (D) €ubic Conve hutaroa 





-2.0 -4.0 0.0 1.0 2.0 


(E) Piecewise Linear (F i Sine@tune ten 


Figure 3.1 Interpolation Functions(1D). 
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sampled data values at the interpolating knots or sampled 
Pouies mei CbNer swords, 91f£ F 1s a2 sampled function and G is 
the corresponding interpolating function, then G(Xk) = F(Xxk) 
whenever XK is an interpolation knot. For equally spaced 
data samples, many interpolation functions can be written in 


the form 


eile) “4% On (see ) (cous) 


where h represents the Sampling increment, the Xx s are the 
interpolation knots,and U is an interpolation kernel. The 
interpolation kernel in Equation 3.1 converts discrete 
data into a continuous function by an operation similiar to 
eme convolution. Mathematically developed interpolation 
kernel functions include nearest neighboring, linear, cubic 
weeame, cuble convolution, piecewise linear, and sinc func- 
on. ame adimensronal (iD) sversion of these interpolation 
mimetrons are illustrated in Figure 3.1. [ese. she The 
mirst funetion in Figure 3.1 1S asample and hold, nearest 
Meaehboring, or replication interpolation. The second func- 
tion iS a Square convolved with another square which results 
mama Criangle or linear interpolation. haem e lal ad etune ton 
is the convolution of three squares or a cubic spline. The 
meuUGeh fLuUmction iS an arbitrary piecewise linear function 
which illuminates the arbitrariness of the interpolation 
meme te tie fitth function is a cubic convolution kernel 
which is composed of piecewise cubic polynomials defined on 
the subintervals (-2,-1), G1 Oi. On le ATi Gagmmiclen 2 ) 
beef. 3]. The last one is the sinc function which provides 
mae Exact reconstruction. The use of the sinc function 
Meamx/X) as am interpolant, which has negative valued side 
lobes, requires an infinite number of terms. Realization is 
ret. cult. eer et (een thie  /red ligation. Polynomials of 
order two or greater can also be employed as an interpola- 


tion kernel function. The amplitude spectra of the nearest- 
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CUBIC CONVOLUTION —=——- 
LINEAR INTERPOLATION -—-— 
NEAREST-NEIGHBOR ———— 


AMPLITUDE 





0 1/h 211/h 31/h 4th 
FREQUENCY 
Pieures3c7 Amplitude Spectra of Interpolation Functions. 





(A)Nearest Neighboring  (B)Linear Interpolataen 


-eso 
7 


Recanstructed signal me 


EM 


——— 
 E_— NO 





CC caries Punic tabon 
Figure 3.3 One-Dimensional Interpolation Process. 
neighboring, linear interpolation, and cubic convo}@iwem 


interpolation kernels are shown in Figure 3.2 for frequency 


from O to 4T/h (where h is the sampling interval 


36 


[Ref. 3]. The response of an ideal interpolation kernel is 


a unit step function inthe frequency domain. In image 
data, the loss of high frequency information causes’ the 
image to appear blurred. On the other hand, deviation from 


the ideal spectrum beyond the shaded area contributes’ to 
pela 1118 - One-dimensional interpolation examples with 
several interpolation kernels are performed in a fashions 
gown in Figure 3.3. [Ref. 4]. Interpolation kernels have 
a Significant impact on the numerical behavior of the inter- 
polated functions. Because of their influence on accuracy 
and efficiency, it iS neceSsary to select carefully an 


appropriate interpolating kernel for an image processing. 


2. Two-Dimensional Interpolation 





(A)Piecewise Linear Interpolation (B)Bilinear Interpolation 


Figure 3.4 Two-Dimensional Linear Interpolation. 


A two-dimensional(2D) interpolation is accomplished 
by two seperate 1D interpolations with respect to each coor- 
dinate. It should be performed along separable orthogonal 
coordinates of the continuous signal. A 2D linear interpo- 
lation kernel function is an example of an orthogonally 
separable interpolation function: G(x,y) = G(x)*G(y). 


Faeure 3.4 shows two examples of the linear interpolation 
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method in 2D sittimeons: bse ce 4) Example (A) shown in 
Figure 3.4 is performed in a piecewise fashion. In region I 
of example (A), points are linearly interpolated in the 
plane defined by pixels A,B,C, while in region II, the 
interpolation is performed in the plane defined by pixels 
Jah e (Ore 1b): The continuous bilinear interpolation is “sh@viaeem 
example (Bae It is done by linearly interpolating points 
along separable orthogonal coordinates of the continuous 


image Signal. 


C. PROBLEMS IN 3-DIMENSIONAIS DATA EEO eaves 


A 3D data is a collection of discrete values resulting 
from equalspace sampling. This data consists of small iden- 
tical parallelepipeds (voxels) divided by three sets of 
planes parallel to the K, Y, and zaxes. Each voxel hase 
Sample value. It is referred to as the intensity value of a 
voxel. For clarity of discussion the voxel value is e@eamee 
on the grid of the object-coordinate system. Figure 3.5(A) 
Shows the relationship between the voxels and the object- 
coordinate system before rotation. Orthogonal projectaan 
along the Z-axis is performed by integrating the intensity 
value of voxels along a line parallel to the Z-axis( any 
point P lying inone of the X-Y plane in the obgie@eme 
coordinate can not be obscured with another point on the 
same plane). The object-coordinate grids exactly correspond 
to the viewing-coordinate grids. 

After rotation, the object-coordinate grid does not share 
the same location as the viewing-coordinate grid as shown in 
Peo ET O2 75 75) (bn) The coordinate value of a voxel on the grid 
of the object-coordinate will be transformed to its Viewines 
coordinate values. These values may not be integers which 
are the viewing-coordinate grids. This is the reason why we 
have to interpolate intermediate intensity values at the 


Viewlne- Coordinate erider 


ers) 






x 
Viewing- 
Coordinate 
Ne 
(A) (B) 
Mibewre 3.) | ee aoe between 3-Dimensional Data and 
Viewing-Coordinate Before Rotation(A) and After Rotation(B). 


To reconstruct the exact object, a high order interpola- 
tion kernel such as sinc functions will produce better 
results. But uSing a high order kernel requires many calcu- 


lations and sometimes it exceeds reasonable computer capa- 


Paslity. Lf we can achieve an appropriate result with a low 
erder kernel function, We can realize a trade-off between 
fidelity and processing time. For that reason, only the 


linear interpolation kernel in the 2D and 3D cases will be 


used in this work. 


D. APPLICATION OF LINEAR INTERPOLATION 
a. One-Dimensional Interpolation 


PEE imMedihaEncetpolation function will yield lower 
quality results compared to that of the high order interpo- 
melon function . An intermediate interpolated point is 
calculated with the nearest two sample points. Linear 
interpolation includes an assumption that the two adjacent 


sample points have a linear relationship. 


Sy, 


Let us consider a linear interpolation examp lemme 
signal after rotation in ae eumeace. Figure 3.6 illuminates 


this Siltuatwonme 


Be 


%A Object-Coordinate 









6 Rotation 


A 
X Viewing-Coordinate 


Figure 3.6 Projection Process after Rotation in 
One-Dimens Pome 


Suppose an arbitrary function is rotated in the Commeems 
clockwise direction by an amount 9 degrees. In this case the 


Sample: “podmese © err on cies) will be projected onto the 


viewing-plane (O-A®*). A projected viewing-coordinate value 
at the grid of viewing-plane will have nonintegers. In the 
above figure, xl,x2,x3,°°* are projected onto viewing- 


coordinate values with the sample values (P1,P2, P3geesueem 
Therefore we have to determine the intensity values at the 
the prids of object -coordinace. Because the adjacent two 
sample values have a linear relationship, we can determine 
the intermediate values from these two adjacent sample 


points (Figure 3.6) 


Ple = Pl * (x2 - xl) + P2 * (xl*]- aes 
P2° = P2 * (x3 - x2°) + P3 * (x2¢ =e?) ieee 
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Po owoOCmeePmeeamat ton 3.2 and 3.3, each interpolation 


Peimireseemmbbrttolwecatrtons, 1 addition, and 2 subtractions. 
2. Two-Dimensional Interpolation 


The two-dimensional interpolation iS an extension of 
a 1D interpolation because the Signal function is defined in 
Bie two orthogonal coordinates as G(x,y) = G(x)G(y). Here 
three algorithms are tested with a sinusoidal 2D function 
having 64? elements: ewe =e eas (x+y) + oh27 . For iis 
example an identical sampling interval in both X, and y 
coordinates are used. The center of rotation is always 
taken to be the center of the coordinate systems. The rota- 
tion is applied to the function which is written in object- 
eoordinate(X,Y). After rotation the sample points of the 
@ereginal funetion mo longer coincide with the egrid of 
viewling-coordinates. It iS mecesSary to interpolate the 


Signal value at the grids of the viewing-coordinate. 
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Figure 3.7 Two-Dimensional Rotation. 
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a. Interpolation by Distributing the Sample Vata 


to the Grids of the Viewing-Coordinate 


Since the data values are samples on the grids 
of the object-coordinate, interpolation of the intermediate 
signal should be done in the object-coordinate. But for the 
practical reason of data storage the interpolation by 
density distribution in the viewing-coordinate is attemptéd 
feta Siar In this case, a linear interpolation kernel on the 
viewing-coordinate is assumed similiar to the kernel on Eine 
object-coordinate. The algorithm for this method is 

(1) Consider one object-coordinate grid. 

(2)  Transtorm “at ema a viewing-coordinate value by 
coordinate rotation. 

(3) Distribute a pixel value inversely proportional toa 
the distance from the nearest 4 Surrounding viewing- 


coordinate grids. (see Figure 3.8(a)) 


Object -Coordinate 


Pixel Valus Assumed Coordinate 


Parallal to Viewing- 
Coordinate 











1(ix,iy) \IPYCixtl iy) 


(A) (B) 


Figure 3.8 Interpolation Kernel over Viewing-Coordinate(A) 
Sampling Value Distribution (ee 
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Pee Om omeeoey<- 1x) )][1-(vyery) |} 
eee | Pome la -ix))|lvy-iy] 
eee Po (yx aee) | [ (vy-iy)] 

Ee =*[ Po ~ (vx-1ix)|[1-(vy-iy)] 


Where Ple,P2e ,P3° ,P4* are distributed components of 
the 4 neighbors of the original sample Po located at 
(vewvy). The displayed image thus obtained is shown in the 
Paeteeesoe9. this fiepure indicates Ehat the error is propor- 
tional to the rotation angle and reaches a maximum at 45 
degrees rotation. As shown in the 1D interpolation case, 
the pixel value of the rotated object-coordinate grid will 
be projected orthogonally onto the viewing-plane (x*-ye 
plane). Even though the object coordinate grids are evenly 
Spaced, the projected grids become compressed on the 
viewing-plane. But this method assumes the evenly spaced 
grids regardless of a rotation angle. Since the amount of 
error introduced depends on the _ rotation angle, ist. 2S 


impractical to apply it in the real display. 


b. Interpolation of the Signal from the Nearest 4 
Samples at the Grids of the Object-Coordinate 


In this algorithm, interpolation is performed in 
bilinear fashion on the object-coordinate. The algorithm 
for this method is 

(1) Consider one viewing-coordinate prid. 

(2 ) Transform it into object-coordinate value by an 
mowerse coordinate rotation. 

(3) Peruonumime thes tneerpolation. -from its nearest 4 
Surrounding samples by the bilinear method. 

This is a correct linear interpolation. However problems 
with this algorithm are (1) We have to rotate all viewing- 


coordinate grids. Ce) All input data should be accessed 
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(CO) 30° Degrees Rotation (D) 45 Degrees Rotation 


e 


by DES aoe e Sample Values to the Grids 
of the Viewing-Coordinate. 


Piel e alse 2 arte: Image with the Interpolation 


Simultaneously in computer memory, which requires’ large 
working space. 
Even with a large computer system, individual working 


Space 1S not usually large enough for storing all data 
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Pixel Value 





Object-Coordinate 


Object -Coord {nate Viewing-Coordinate 


(A) (B) 
Figure 3.10 Interpolation Kernel over te ee (A ) 
anemepitinear Interpolation(B). 


values. Therefore this method is. not Pprueeumeas for 
mamepimlation of 3D data. 

| A displayed image obtained by this method is 
Siem in Figure 3.11. It shows a better result than the 
Mime polation after rotation of object-coordinate. Blurring 


effect is due to inherent linear interpolation error. 


Gx Interpolation Using Cone-Shape Kernel 


We have already disscussed two interpolation 
methods. But these two methods have flaws which make it 
peer icult to apply it to the 3D situations. ene Git ed 
method which avoids these flaws uses a 2D cone-shape kernel. 
Figure 3.12 illustrates how the cone-shape kernel is used 
meee D linear interpolation. The density distribution of a 
pixel in the exact 2D linear interpolation should be done 
like Figure 3.12(B). But the cone-shape density distribution 
method is used for this purpose. The advantage of using a 


cone-shape kernel is 
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(6) "30 Derrees Rotat ron (D) 45 Degrees Rotation 


[Qo Uedbinay-wee sans ad Displayed es by eg ae ee nS the Signal 
from the Nearest 4 Samples at the Grids 
of the Object-Coordinate. 


(1) It is invariant to the direction) of "rotagucue 


(2) Interpolation is performed on the object-coordinate. 


“The algorithm for this methods 


(1) Consider one object-coordinate jomiee 
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Pixel Value Pixel Value 





Object-Coordinate 


(A) (B ) 


Bee eo 3. 12 (A) Cone-Shape Kernel (B) Rectangular kernel. 


fetcansform it to the viewing-coordinate value by 
Goomarmate rotation. 

mo Calculate the distance from it to the 4 surrounding 
apyeet-coordinate grids. 

fee it a distance is smaller than unity, calculate the 
distribution component of a intensity value inversely 
proportional to the distance. Otherwise asSign zero. 

eeeoumesthe —distributions of alt poimes to pet the 
intensity value at the grid of the object coordinate. 

Using this kernel, only viewing-coordinate grids inside 
the kernel are affected by the sample value at the center of 
the cone. Otherwise it is ignored. BPaewreno. 13 indicates 


that the error introduced by using the cone-shape kernel is 


mot severe. Resulting displayed image is almost the same as 
the above one. Because this cone-shape kernel is invariant 
memene direction of rotation, it can be extended easily 


into a sphere-shape kernel in a 3D situation. 
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C) 30 Degrees Rotation (D) 45 Degrees Rotation 


Pa pure 63.45 Displayed Images by Interpolation 
Using Cone-Shape Kernel. 


The disadvantages of this method are: (1) ee 
involves many multiplications to calculate the distam@eus 
(2) Error includes inherent linear interpolation error plus 
the cone-shape kernel error. The advantage of this method 
LSe tliat ane ean eaSily be applied in a 3-dimensiomag 


Siruac lone 
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3. Three-Dimensional Interpolation 


NomemeloceGmetinethne 82 Situation, linear interpola- 
tion using a cone-shape kernel is the most appropriate 
method. Because the shape of this kernel iS invariant with 
respect to the direction, it can be easily extended to the 
Speimrerpolation. Applying the linear interpolation method 
to the 3D situation requires avery complicate procedure. 
ereer rotation, a parallelepiped looks like ae polygon 
instead of a cube to the sight of the viewer. Therefore, 
performing linear interpolation to the sphere is easier than 
meine polygon. Lin these -case., only grids of object- 
coordinate inside of the sphere are affected by a pixel 
value lying at the center of a sphere. weal eerathm £or 
this method is 

(1) Consider one object-coordinate grid. 

fom itanstorm 1t into viewing-coordinate value by 
Se9erdainate rotation. 

(3) Calculate distance from it to the surrounding 8 
object-coordinate grids. 

(4) Calculate the distribution of intensity value at the 
8 object-coordinate grids when a distance is smaller than 
unity. Otherwise ignore it. 

(5) Sum the distribution components of the 8 surrounding 
intensity values’ to calculate the value at a object- 
Seecdinate grid. 

The disadvantages of this method are (1) Necessity for 
mime hoagtelae- point miltiplications, (2) Total error is the 
sum of inherent linear interpolation error and sphere-shape 
Memmel error. More details will be discussed in the Chapter 
ae 
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IV. PROGRAM IMPLEMENTATION OF REPROJECTION AND Disse@itae 


Asy PORTHOGRAPHICG REPROIEC Cl Te 


Orthographic reprojection is performed numerically in the 
computer by summing the value of voxels along parallel paths 
through the reconstructed volume. The reprojected image 
constitutes a 2-dimensional image on the screen. 

Figure4-l illustrates the orthographic reprojection 


process onto the X*-Y* plane along Z* coordinate. 











5G 
caer ea 






Summation Paths 


jected 
Reprojecte (Yk th Slice veg 


Pixels on the 


4 
| 
Viewing-Plane |__|! _ Volume ) 
(Yk th Row? Sct 
Figure 4s | Reprojection Process before Rotabaeme 
To implement this process in the program, firs tae 
viewing-plane (2-dimenstomalwancay: "Viewpln') is defined 


parallel to the viewing-coordinate X*-Y* shown as in program 
‘Rot 3 Dim’. The viewing-plane composes of X*(column) and 
Ye(row) coordinates and lies orthogonally to Z*-coordinate 
(X*,Y*,Ze are viewing-coordinates). The Yx'th slice of the 
volume is projected onto the Yx ‘th row of pixels in the 


viewing-plane as shown in Figure 4.1. 
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Before rotation, every voxel is arranged along viewing- 
coordinate soa object-coordinate grids COmmeSpDOndesrto 
viewlng-coordinate grids. 

The reprojection can be performed as in the following 
procedure: 

(1) For every incoming voxel, read the viewing- 
coordinate Ni ehereeatenesum the demnsaty vale onto a pixel 
Smeresponding at XK°,Y¥K° in the viewing-plane (2 Sea ieeay 
Sasewpln ). | 

Viewpln(X«*,Yx*) := Viewpln(Xx*,Yu°) + Pixel (Xk*,Y«*,Z Ko); 

(2) Continue to reach the last voxel. 

After rotation, it is not clear what portion of a voxel 
will contribute the pixel value in the viewing-plane because 
each voxel looks like a polygon to the sight of view. 
Through the interpolation method discussed in the Chapter 
IIL, the appropriate portions of the voxel value are added 
to the pixel values in the viewing-plane (see subroutine 


‘Projection’ and ‘Dissection’ of the program 'Rot_ 3 Dim’). 


Viewing-Plane 








Summation Paths 


Repro jected 
Pixels on thé 
Viewing-Plang 


4 
, | 


Baieure 4.2 Reprojection Process after Rotation. 


jek 


The resulting pixel values in the viewing-plane will 
exceed the allowed pixel values of the display device so 
that the pixel should be rescaled to fit the specific 
display device requirement. 

Several reprojected images with sphere-shape kernel are 
illustrated in the Figure 4.3. Different viewing-angles 
provide more informations about the structurpern of 


3-dimensional data. 


De 








30 Degrees 45 Degrees 75 Degrees 
Layia Iepeshealcjo, clevenive @u: eat 














15 Degrees 90 Degrees 180 Degrees 
(CE Rotation. about Y>Axis 





30-30 Degrees 45-45 Degrees 60-60 Degrees 
Ceomeroratltomeapour boLM-Ax1s 


Figure 4.3 Displayed Images by Reprojection. 
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B. SINGLE PLANE DissreriG 


The dissection is the technique that an object plane is 
clearly displayed showing spatial relationship with other 
planes in the volume data. But only single plane dissection 
1s considered in this work because the dissection pr@eeeee. 
very similar to the single plane dissection. Single plane 
dissection is performed by removing all other planes’ and 
displaying an object plane. 

Suppose we want to see a plane to the selected viewing- 


angle after rotation. An viewing-plane will be lying in the 


—, i 





Figure 4.4 . Single Plane Dissection. 
Viewing-coordinate shown as in Figure 4.4. After Protairione 
vilewlng-plane cuts the volume with an arbitrary angle. The 


voxels reprojected onto the viewing-plane are not squares 
rather polygons with an arbitrary shape. Therefore it is 
difficult to determine what portion of sample values will be 
contributed to the pixels in the viewing-plane. 

To calculate the intensity value of pixels meee 
viewing-plane, 3-dimansional interpolation is necessary. The 


interpolation using sphere-shape kernel is used to calculate 
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the intensity values on the viewing-coordinate grids. 
Subroutine '‘Dissection" in the program 'Rot_3 Dim' shows the 
dissection process. Procedure for this scheme is 

ae Accept a desired rotation angle and desired 
Z*-coordinate to see. 

(2) Take a voxel from the data. 

) Tramotomnm LC el pete the corresponding viewing- 
coordinate value. 

(4) Devide the voxel value to the Surrounding 8 
coordinate-grids by 3-dimensional interpolation method. 

(se it the Z*-coordimate of Ehe grid corresponds to the 
desired Z*-coordinate, add the portion of the voxel value to 
the corresponding pixel in the viewing-plane. Otherwise 
memore it. 

(6) Continue to the last voxel. 

Several dissected planes from the reprojected image are 
miiustrated in Figure 4.5. We can see the checker-board 
effect in some displayed images which does not appear in the 
figure. This means that the linear interpolation with 
sphere-shape kernel is _not sufficient for the single plane 
eerseceeton in the certain rotation angle. If one want 
higher Quality images, more precise and efficient 


interpolation method should be developed. 
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45 (ZCRD=25) 60° (ZCRD=-25 


2 90° (ZCRD=0) 
(A) ROtAErOn abour = xis 





45” (A2GRDe2 On. 135. (GER D= 2m) es el Goi 0) 
(B) Rotation about’ Y-Axis 


45° -45°(ZCRD=0) 
Both-Axis 





Figure 4.5 Deane Image by Dissection 
( s the Z th Slice of Volume Image). 
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VeRPcen@EustON AND RECOMMENDATIONS 


A. CONCLUSION 


An inherent disadvantage of the reprojection method is the 
Sescuring of underlying and overlying structures in the 
2-dimensional viewing-screen. Therefore it is desirable to 
have better visibility.of selected portions of the volume 
data. Interactive user selection of an adequate viewing- 
angle and the application of the dissection method is 
proposed to resolve this problem. Another important 
problem is’ the need to manipulate a huge amount of data 
as efficiently as possible. 

Three different algorithms were studied here for the 
Seordinate transformation: Dire 6 tates 1 1X sei t 1 plea Eon ; 
Decomposition with fixed coordinate, and Recalculation and 
Indexing. 

In the Precalculation and Indexing method, the compo- 
nents of coordinate and matrix element are calculated itera- 
tively and stored in an array. Then these S OMERIRIEIES are 
indexed according to the incoming coordinates and” used in 
maemcalculation so that the number of matrix multiplication 
1s decreased. This method requires an appropriate size of 
memory. It is possible to rotate 64? data elements in 25 
seconds with the VAX-11/750 minicomputer. 

To realize this algorithm in hardware, two hardware 
design schemes have been suggested in Chapter II. 

After rotation to a selected direction, the positions of 
the voxels do not correspond to the viewing-coordinate 
grids. Therefore several linear interpolation algorithms are 
Studied. 

The reason why linear interpolation using a cone-shape 
Memme: Cam be extended from 2D to 3D is attributed to its 


invariant shape with respect to the different directions. 


= / 


The sphere-shape kernel is a 3D extension of the cone-shape 
kernel in the 2D situatience r 

Using this 3D interpolation method, a reprojection and 
Single plane dissection program has been tested with an 
areca DS -saatac Although linear interpolation Wsamme 
sphere-shape kernel has inherent errornse this method 
provides the reasonable quality images in a reasonable 
amount of time. 

Manipulation of a huge amount of data such asia 
computed tomography data requires a continuous’ trade-off 
between processing time and memory, or between processing 
time and image quality. The only way to resolve this imme 
of problems is to use heuristic appuetenoce First, select 
an initital algorithm and test it. If the result "seine 
adequate for the specific purpose, try another algorithm. 
This procedure is continued until an appropriate result can 


be reached. 


B. RECOMMENDATIONS 


The coordinate rotation and reprojectiton method can be 
widely used in the viewing or manipulating of 3-dimensional 
data. The algorithms which have been developed in this work 
may not be optimum because these are the result of heuristic 
approaches. 

The coordinate transformation algorithm implemented here 
can rotate any reasonable size of 3D data in a short time 
But, further studies are required in order to get better 
quality images from a 3D interpolation algorithm. 

It is worth while to test the implementation of the 
reprojection algorithm with real CT scan data. in EhiSgueee 
new problems may arise, which will require further studies. 

The implementation of the dissection and the dissolution 
capabilites are very helpful for the investigation of 


Structures am 3D volume meaece 
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